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Abstract 

The inclusive photon energy spectra measured by the Large Hadron Col- 
lider forward (LHCf ) experiment in the very forward region of LHC proton- 
proton collisions at \/s = 900 GeV are reported. The results from the analysis 
of 0.30 nb^ 1 of data collected in May 2010 in the two pseudorapidity regions 
of rj > 10.15 and 8.77 < rj < 9.46 are compared with the predictions of 
the hadronic interaction models DPMJET 3.04, EPOS 1.99, PYTHIA 8.145, 
QGSJET 11-03 and SIBYLL 2.1, which are widely used in ultra-high-energy 
cosmic-ray experiments. EPOS 1.99 and SYBILL 2.1 show a reasonable 
agreement with the spectral shape of the experimental data, whereas they 
predict lower cross-sections than the data. The other models, DPMJET 3.04, 
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QGSJET 11-03 and PYTHIA 8.145, are in good agreement with the data be- 
low 300 GeV but predict harder energy spectra than the data above 300 GeV. 
The results of these comparisons exhibited features similar to those for the 
previously reported data for yfs = 7TeV collisions. 
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1. Introduction 

The observations of ultra-high-energy cosmic rays (TJHECR) have made 
notable improvements in the last few years [H Q j^] Q p] 0] Q ■ However, al- 
though some critical parts of the interpretation rely on the Monte Carlo (MC) 
simulations of the air shower development, very forward particle emission in 
the hadronic interactions, which are relevant to the precise understanding of 
air showers, have been poorly understood thus far, especially at such high 
energies. To reduce the uncertainty in MC air shower simulations, the Large 
Hadron Collider forward (LHCf ) experiment has performed measurements of 
the neutral particles emitted to the very forward region of proton-proton col- 
lisions at the LHC. In 2010, the operations at a/s = 7TeV and 900 GeV were 
completed. The photon energy spectra obtained from the 7TeV data have 
been previously reported elsewhere [1] . Below this energy, the measurement 
of the Pt spectra of 7r°s by UA7 in the rapidity region Y = 5.05 - 6.65 is 
available. It is interesting to have a single-photon measurement at degrees 
at different collision energies in order to discuss energy dependence. In this 
paper, we report the inclusive photon energy spectra in the very forward re- 
gion for ^/i = 900 GeV proton-proton collisions with the same detectors and 
analysis methods as those used for the y/s = 7 TeV analysis. 

Two LHCf detectors, called Arml and Arm2, were installed in the instru- 
mentation slots of the TANs (Target Neutral Absorbers) located at ±140 m 
from the ATLAS interaction point (IP1) and covering the pseudorapidity 
range from 8.7 to infinity (zero degrees). Each detector had two sampling 
and imaging calorimeters composed of 44 radiation lengths (1.55 hadron in- 
teraction lengths) of tungsten and 16 sampling layers of 3 mm thick plastic 
scintillators. The transverse sizes of the calorimeters were 20 mm x 20 mm 
and 40 mm x 40 mm in Arml and 25 mm x 25 mm and 32 mm x 32 mm in 
Arm2. The smaller and larger calorimeter of each Arm are called the small 
tower and large tower, respectively. The cross sections of the calorimeters, 
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as observed from IP1, are illustrated in Fig.[TJ During the operations that 
were used in the analysis reported in this paper, a large fraction of the large 
Arml tower was obscured from the IP as indicated by the shaded area in the 
Fig.[T]due to the beam pipe material between the IP and the detector. 

Four X-Y layers of position-sensitive detectors (scintillating fiber, SciFi, 
belts in Arml and silicon micro-strip sensors in Arm2 with 1-mm and 0.16- 
mm readout pitches, respectively) were inserted to measure the transverse 
positions of the showers. The LHCf detectors have energy and position 
resolutions better than 5% and 200 /im, respectively, for >100 GeV pho 



tons. Detailed descriptions of the detectors can be found elsewhere [10f| [11 
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This paper describes the first results of the analysis of inclusive photon 
energy spectra for -y/i = 900 GeV proton-proton collisions, which are primar- 
ily produced from the decay of 7r° and rj mesons generated in the collisions. 
The data set and the MC simulation used in the analysis are introduced in 
Sec. [2] and Sec.[3j respectively. The analysis process is described in Sec.HJ 
The experimental results and comparison with the MC predictions of several 
hadronic interaction models are presented in Sec.|5]and summarized in Sec. El 

2. Data 

The data sets used in the analysis were taken on 2, 3 and 27 May 2010 
during the LHC operations with proton-proton collisions at a/s = 900 GeV, 
which correspond to the LHC fill identification numbers (Fill ID) 1068, 1069 
and 1128, respectively. In these fills, the LHC operated with one crossing 
bunch and one non-crossing bunch at IP1 in Fill IDs 1068 and 1069 and 
with four crossing bunches and three non-crossing bunches at IP1 in Fill 
ID 1128. The luminosity (L) at IP1 during these fills was measured by 
the ATLAS experiment [16|. The luminosity during Fill ID 1068 and 1069 
were 8 - 3 x 10 27 cm~ 2 s~ 1 and 12 - 4 x 10 27 cm _2 s _1 , respectively. The total 
luminosity of the four crossing bunches in Fill ID 1128 was approximately 
8 x 10 27 cm _2 s _1 . The total integrated luminosity (J Ldt) during the LHCf 
operations in the three fills was 0.30 nb _1 . The uncertainty of the lumi- 
nosity determination is ±21% [16]. The inelastic cross-section (ai ne i) for a 
y/s =900 GeV proton-proton collision was estimated to be <7i ne i =53.0mb 
from the predictions of the total cross-section and the elastic cross-section, 



which are based on the recent experimental results [17J |18[ . The number of 
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Date 


Time (UT) 


Fill ID 


L (cm 2 s 1 ) 


j Ldt 


Crossing bunch 


2 May 2012 


12:50 - 19:23 


1068 


8~3x 10 2Y 


0.11 nb^ 1 


lxl 


3 May 2012 


00:17 - 07:08 


1069 


12~4x 10 27 


0.17 rib" 1 


lxl 


27 May 2012 


13:18 - 14:03 


1128 


8 x 10 27 


0.02 rib" 1 


4x4 



Table 1 : Summary of the luminosity during y/s — 900 GeV operations [l6j . 



inelastic collisions (iVj ne during the three fills was calculated to be 1.58 x 10 7 . 
The luminosity during the three fills is summarized in Tab {I] 

During the \fs = 900 GeV fills, the LHCf operations were performed with 
a high-gain operation of the PMTs for the sampling layers to detect photons 
with energies as low as 50 GeV with a nearly 100 % trigger efficiency and 
with a lower threshold with respect to the \/s = 7 TeV data. The typical 
PMT gain for the high-gain operations was 3-5 times higher than the nom- 
inal gain that was used to obtain the 7 TeV data. Neither the saturation 
of PMTs nor the range of ADC caused problems because the maximum en- 
ergy of the incident photons were expected to be 8 times lower than those 
for the 7 TeV data. The data acquisition (DAQ) triggers were generated 
from beam pickup signals (BPTX) followed by "shower trigger" signals. The 
trigger condition of the shower trigger was that signals from any three suc- 
cessive scintillator layers in any calorimeters exceed the predefined threshold 
(approximately 17 MeV for high gain). The average DAQ live times during 
the LHCf x/s =900 GeV operations was 99.2% (Arml) and 98.0% (Arm2). 
The total numbers of triggered events in Arml and Arm2 were 44,389 and 
62,916, respectively. Because of the very low luminosity and the low event 
rate per inelastic collision, the probability of the pile-up of events was < 10~ 4 , 
negligibly small. 



3. MC simulation 



To compare the experimental results with the predictions of hadronic 
interaction models, MC simulations were performed with the hadronic inter- 

SIBYLL 2.1 



action models, QGSJET 11-03 20[, PYTHIA 8.145 |21| |22j 



EPOS 1.99 241 and DPMJET 3.04 25 . In the MC simulations, 3 x 10 7 
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inelastic proton-proton collisions were generated by each model, and the sec- 
ondaries were transported in the beam pipe from IP1 to the LHCf detectors. 
The magnetic fields of the dipole magnets located between IP1 and the LHCf 
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detectors were taken into account. The detector response was calculated us- 



ing the EPICS 8.81/COSMOS 7.49 simulation package [19 . 

In addition, 1.0 x 10 8 events were generated using QGSJET 11-03. This 
data set was used for studies of the detector response and particle identifica- 
tion (PID) correction described in Sec. 14.21 

4. Analysis 

4-1. Energy Reconstruction 

The sum of the energy deposited in the 2 nd to 13 th scintillator layers, after 
corrections for gain variation and the non-uniformity of the light yield of each 
scintillator layer, was used as an energy estimator for the primary photons 
incident on the LHCf detectors. Each PMT gain was premeasured for various 
HVs using N2 laser calibration. The calibration of the deposited energy in 
each layer was performed for different HV settings (low-, nominal-, and high- 
gain operations) using 50 - 200 GeV/c electron beams and 150 GeV/c muon 



beams at the CERN SPS [131 ] . The non- uniformity of the light yield of each 
scintillator layer was measured using a /3-ray source before assembling the de- 
tectors. Because a fraction of the shower particles leak out of the sides of the 
calorimeters ('shower leakage'), the total energy deposited was corrected for 
'shower leakage' by a function of the shower impact position. This function 
was determined by the MC simulation. The impact positions of the show- 
ers were determined using the information from the position-sensitive layers. 
The events that fell within 2 mm of the edges of calorimeters were removed 
from the analysis to avoid the degradation of the energy resolution due to 
'shower leakage'. We set the energy threshold of this analysis to 50GeV to 
avoid background from the interactions between secondary particles and the 
beam pipe, which was expected to be concentrated below 50 GeV according 
to a MC simulation. The trigger efficiencies of both Arms were also checked 
by two samples. One was an unbiased data sample with nominal gain trig- 
gered by the shower triggers at the opposite side of the detector. The other 
was the detector MC simulation for both high and nominal gains. Consid- 
ering the difference of gains (a factor 3-5) between the high and nominal 
gain operations, the two methods gave consistent results. We found that the 
efficiencies were 100% for > 30 GeV incident photons for both arms for the 
higher gain. This was sufficiently lower than the 50 GeV analysis threshold. 

Similar to the previous 7 TeV analysis, the systematic uncertainties of the 
absolute energy scale were evaluated from the reconstructed invariant mass 
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Energy range (GeV) 


50 - 150 


150 - 200 


200 - 250 


250 - 300 


300 - 450 


Arml Small (%) 


-2, +3 


-6, +3 


-26, +5 


-48, +14 


-48, +14 


Arml Large (%) 


-8, +5 


-20, +5 


-20, +5 


-50, +15 


-71, +30 



Table 2: Summary of the systematic errors for the energy spectra of Arml due to the 
energy scale uncertainty. 



Energy range (GeV) 


50 - 150 


150 - 225 


225 - 300 


300 - 450 


Arm2 Small (%) 


-15, +8 


-15, +8 


-28, +12 


-66, +34 


Arm2 Large (%) 


-5, +7 


-20, +7 


-24, +10 


-63, +14 



Table 3: Summary of the systematic errors for the energy spectra of Arm2 due to the 
energy scale uncertainty. 



of the 7r°s from the 7TeV data taken with a nominal gain. Additionally, 
we checked 7r° mass peaks in the 7TeV data with a higher gain taken in a 
different period. We found +2.7% (+0.7%) differences of the mass peaks 
from those for the nominal gain in Arml (Arm2). These differences were 
compatible with the uncorrelated energy scale errors (±3.5 %) quoted for 
the energy scale calibration of the detectors using the SPS beams or a long- 
term time variation. Conservatively, they were added to the energy scale's 
systematic error in quadrature. Finally, the energy scale uncertainties were 
estimated to be [-10.2 %, +1.8 %] and [-6.6 %, +2.2 %] for Arml and Arm2, 
respectively. The systematic uncertainties of the energy spectra due to the 
energy scale uncertainty are listed for the Arml and the Arm2 detectors in 
Tab. [2] and Tab. El respectively. 

The energy of photons in multi-hit events with more than one photons 
incident on a single tower would not be reconstructed correctly. We estimated 
the possible bias in the energy reconstruction due to double-incident events 
using MC simulation using QGSJET 11-03. We found that the number of 
events whose true energies were modified by more than 2 % was very small. 
The fraction of such events was expected to be less than 1 % of the events 
having a total incident energy of >40 GeV. Although such multi-hit events 
can be identified using the lateral distributions measured by the position- 
sensitive layers, we did not apply the multi-hit cut or any correction to the 
spectra in this analysis to avoid a bias due to the misidentification of single 
photon events as multi-hit events. The effect was at most 1% in the lower 
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energy bins. 

4-2. Photon Event Selection 

To select the electromagnetic shower events and to eliminate contamina- 
tion by hadronic shower events, a parameter called L 90 % was defined. L 90 % 
is a longitudinal length in units of radiation length (r.l.) in which 90 % of the 
total shower energy is deposited in the calorimeter. Figure[2] shows the L 90 % 
distribution of the Arml small tower events with reconstructed energy in the 
range of 50 to 100 GeV. The two peaks near 13 r.l. and 35 r.l. correspond to 
electromagnetic showers and hadronic showers, respectively. Fig. [2] also shows 
the L 90 % distributions generated by the MC simulation with QGSJET 11-03 
for pure photons and pure hadrons (neutrons). These MC distributions have 
been normalized to the L 9 o% distribution of the experimental data. They are 
hereafter called the 'template'. We set the L 9Q % criteria to keep the pho- 
ton selection efficiency epiD = 90% over the entire energy range based on the 
template for photons. In Fig.[2]this would correspond to L 90 % = 16.8 r.l. 

The purity (V) of a photon sample was estimated by normalizing the 
templates for photons and for hadrons to the measured L 9 o% distribution for 
each energy range ('template fitting' of the L 90 % distribution). V was defined 
as V = N photon / (Nphoton + Nhadron) in each energy region. Here N photon and 
Nhadron are the numbers of photons and hadrons in the selected L 99 % range 
in the template, respectively. The correction factor V x e P ] D was applied to 
the number of events in each energy bin to correct for the inefficiency of the 
photon selection and for the residual contamination by hadrons. 

However, there were small discrepancies between the L 90 % distributions of 
the experimental data and of the MC simulations. These discrepancies may 
be caused by errors in the absolute energy determination or in the channel- 
to-channel gain calibrations. Here, we consider the systematic uncertainty 
caused by the uncertainty of the 'template fitting' method for obtaining pho- 
ton spectra. Small modifications of the template (widening with respect to 
the peak position up to 30 % and a constant shift up to 1.0 r.l. for Arml and 
0.8 r.l. for Arm2, to give the best match with the data) provide the size of 
uncertainty in the correction factors to the photon spectra. The difference of 
the correction factors between the original and the modified template meth- 
ods amounted to ±10% (±12%) and ±55% (±44%) for below and above 
150 GeV photon energy, respectively, in the Arml small (large) tower. The 
difference of the correction factors were ±20 % (±25 %) in the Arm2 small 
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(large) tower for entire photon energy. These numbers were assigned as the 
systematic uncertainty of the energy spectra due to PID errors. 



4-3. Background Subtraction 

The background particles from the interactions between the proton beams 
and the residual gas in the vacuum beam pipe hit the detectors synchronously 
with the beam-beam events. The amount of this background can be esti- 
mated using the events triggered by the passage of non-crossing bunches. 
Assuming that the beam intensity of each bunch was same, the background 
levels were estimated as approximately 1 % and 2 % for the small and large 
towers, respectively, in both Arms. The estimated backgrounds were sub- 
tracted from the energy spectra. 

4-4- Beam Center Position 

The projected position of the zero-degree collision angle at the LHCf de- 
tectors, referred as the 'beam center', is an important parameter in the geo- 
metrical analysis of the experimental data. Because the flux of the secondary 
particles produced by = 900 GeV proton-proton collisions was expected to 
be uniform over the acceptance of the LHCf detectors, the beam center could 
not be determined directly from our measurements. In this analysis, we as- 
sumed that the beam center was at the center position determined by the 
alignment survey of the detectors. The beam center was located near the 
center of the small calorimeter of each Arm as shown in Fig.[TJ The beam 
position and the beam angle were monitored by the Beam Position Moni- 



tor (BPMSW) installed ±21m from IP1 26]. The fill-by-fill fluctuation of 
the calculated 'beam center' at the LHCf detectors during the period for 
the presented data set was approximately 4 mm. We assigned a systematic 
uncertainty of ±2 mm to the beam center. To estimate the effect of this 
uncertainty on the energy spectra, we defined an area that had a slightly 
narrower acceptance than the calorimeter, and an energy spectrum was gen- 
erated from the events falling within this area. Additionally, four spectra 
were made by shifting the area by ±2 mm vertically and horizontally. As a 
systematic uncertainty of energy spectra due to the uncertainty of the beam 
center, we assigned the differences of the spectra as shown in Tab.H] and 
Tab. [5] for the Arml and the Arm2 detectors, respectively. 
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Energy range (GeV) 


50 - 150 


150 - 250 


250 - 450 


Arml Small (%) 


-7, +2 


-7, +2 


-17, +2 


Arml Large (%) 


-11, +6 


-11, +12 


-15, +12 



Table 4: Summary of the systematic errors for the energy spectra of Arml due to the 
beam center uncertainty. 



Energy range (GeV) 


50 - 175 


175 - 250 


250 - 450 


Arm2 Small (%) 


-3, +3 


-4, +9 


-4, +13 


Arm2 Large (%) 


-3, +4 


-3, +6 


-16, +14 



Table 5: Summary of the systematic errors for the energy spectra of Arm2 due to the 
beam center uncertainty. 



5. Energy Spectra Results 

5.1. Reconstruction of Energy Spectra 

To reduce a possible pseudorapidity (rj) dependence when comparing and 
combining the energy spectra measured by the two Arms, we selected Arm2 
events with a pseudorapidity range similar to that of Arml. For the small 
tower, we selected events with the distance (r) from the beam center less 
than 11mm, which corresponded to the pseudorapidity range of rj > 10.15 
(the circles in Fig.[T|). Similarly, for the large tower, we set the conditions 
as 22 mm < r < 44 mm, which corresponded to the pseudorapidity range of 
8.77 < f] < 9.46 (the arcs in Fig. [I]). The calorimeters did not uniformly 
cover the pseudorapidity ranges as shown in Fig.[TJ We confirmed that there 
was a negligible pseudorapidity dependence of the energy spectra inside each 
pseudorapidity range. The reconstructed photon energy spectra of Arml 
and Arm2 are shown in Fig.|3]in units of differential cross-sections da/dEdQ, 
where E is the photon energy and Q is the solid angle. The differential cross- 
section was calculated as da/dEdQ = 1/(J Ldi) dN/dEd(cos9)d(fi, where N 
is the number of events in each energy bin, J Ldt is the integrated luminosity 
after the correction for the DAQ live time for each Arm and 9 and <fi are the 
polar and the azimuthal angles with respect to the beam axis, respectively. 
Considering the geometrical acceptance of the calorimeters, the averages of 
the polar angle (8) are 39 /zrad and 234 /irad for the small and the large tower, 
respectively. The error bars in Fig. [3] indicate the statistical uncertainty; the 
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hatched areas show the systematic errors in which the particle identification 
and the beam position uncertainties were taken into account. Because the 
systematic uncertainties due to the energy determination may be correlated 
between Arml and Arm2[8|, they were not taken into account in Fig.[3j 
The two spectra from Arml and Arm2 in each pseudorapidity region gave 
consistent results within the statistical and the systematic errors. 

The combined energy spectra of Arml and Arm2 are shown in Fig.H] as 
weighted averages, with the weights taken to be the square of the inverse 
of the errors in each energy bin. The error bars of the data (black points) 
represent the statistical error; the hatches in the spectra represent the total 
uncertainty (quadratical summation of the statistical and the systematic er- 
rors). The sources of the systematic error are the particle identification and 
the beam position uncertainties. The energy scale errors were also included, 
assuming a correlation between the two Arms. Note that the uncertainty of 
the luminosity determination (±21 %) is not shown in Fig.HJ It can introduce 
a constant vertical shift of the spectra, but it cannot change the shapes of the 
spectra. We see a smooth spectrum from each of the two pseudorapidity re- 
gions, considering the errors. The similarity of the two spectra suggests only 
a small pseudorapidity dependence between the two pseudorapidity regions. 

5.2. Comparison with Models 

In Fig.HJ the predictions of the hadronic interaction models, QGSJET I 
1-03, PYTHIA 8.145, SIBYLL 2.1, EPOS 1.99 and DPMJET 3.04, are also 
shown. The same analysis processes were applied to the MC simulations as 
to the experimental data except for the particle identification using L 90 % and 
its correction. For the analysis of the MC simulations, the known particle 
type was used. For better visibility, only the statistical errors for DPMJET 
3.04 (red points) are shown by the error bars. Figure^] shows the ratios of 
the MC spectra divided by the data in each energy bin. In Fig.|5j the statis- 
tical error of each MC is shown as the error bar of each point. The trends of 
the experimental data compared to each MC are similar for the two pseudo- 
rapidity ranges. EPOS 1.99 and SYBILL 2.1 show a reasonable agreement 
with the spectral shape of the experimental data, whereas they predict lower 
cross-sections than the data. The other models, DPMJET 3.04, QGSJET I 
1-03 and PYTHIA 8.145, are in good agreement with the data below 300 GeV 
but predict harder energy spectra than the data above 300 GeV. The trends 
of the experimental data compared to the MC predictions in Fig. are similar 
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to those for the single-photon energy spectra in the pseudorapidity r] > 10.94 
previously reported for y/s = 7TeV proton-proton collisions s(]. 



6. Summary 

LHCf measured the forward inclusive photon energy spectra for y/s = 
900 GeV proton-proton collisions in May 2010. The total integrated luminos- 
ity of the data set used in this analysis is 0.30 nb _1 . The two LHCf detectors 
(Arml and Arm2) gave consistent results within the statistical and system- 
atic errors for the small and the large towers, which cover the pseudorapidity 
ranges of rj > 10.15 and 8.77 < rj < 9.46, respectively. The combined en- 
ergy spectra of Arml and Arm2 were compared with the predictions of five 
hadronic interaction models, DPMJET 3.04, EPOS 1.99, PYTHIA 8.145, 
QGSJET H-03 and SIBYLL 2.1. EPOS 1.99 and SIBYLL 2.1 reproduce 
well the shape of the experimental energy spectra, but they predict a lower 
cross-section than the LHCf data. The other models predict harder spectra 
than the LHCf data above 300 GeV. These results of comparison exhibited 
features similar to those for the previously reported data for y/s = 7TeV 
collisions. 
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Figure 1: The cross-sections of the calorimeters viewed from IP1, left for Arml and 
right for Arm2. The cross marks on the small calorimeters indicate the projections of the 
zero-degree collision angle onto the detectors ('beam center'). The shaded areas in the 
upper parts of the figure indicate the shadows of the beam pipes located between IP1 and 
the detectors, where the detectors are insensitive to the detection of IP1 proton-proton 
collision products. The dashed squares indicate the border of a 2 mm edge cut described in 
Sec. 14. II The circles and the arcs indicate the distance (r) from the beam center of 11 mm, 
22 mm and 44 mm, and the pseudorapidity (rj) of 10.15, 9.46 and 8.77, respectively. In 
this analysis, the events in the regions of r < 11mm and 22mm< r <44mm were used. 
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5 10 15 20 25 30 35 40 

Radiation Length[r.l.] 



Figure 2: The Lg % distribution measured by the Arml small tower for the reconstructed 
energy range of 50 - 100 GeV. The black points show the experimental data; the red 
and blue histograms are the L 90 % distributions of MC calculation for pure photons and 
pure hadrons ('templates'), respectively. These template histograms are independently 
normalized to give the best match to the experimental result. 
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Figure 3: Photon energy spectra measured by Arml (red circles) and Arm2 (blue rect- 
angles) in units of differential cross-section da/dEdil. The left and the right panels show 
the results of the small towers and the large towers, respectively. The pseudorapidity cov- 
erages of the small and the large towers are r\ > 10.15 and 8.77 < rj < 9.46, respectively. 
Considering the geometrical shape of the calorimeters, the averages of the polar angle (9) 
with respect to the beam axis are 39 /xrad and 234 /irad for the small and the large towers, 
respectively. The error bars indicate the statistical errors, and the hatched areas indicate 
the total systematic uncertainties that come from particle identification and the 'beam 
center' position. 
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Figure 4: Combined Arml and Arm2 photon energy spectra compared with MC pre- 
dictions. The data from Arml and Arm2 correspond to the integral luminosities of 0.30 
and 0.29 nb _1 , respectively. The left and the right panels are the results of the small 
(i] > 10.15) and the large (8.77 < r/ < 9.46) towers, respectively. The black points 
indicate the experimental data with the statistical uncertainty (error bars) and the to- 
tal uncertainty, quadratical summation of the statistical and the systematic errors (black 
hatches). The systematic uncertainty of the luminosity determination (±21%) is not 
taken into account in the errors. The colored points indicate the results of MC predic- 
tions, QGSJET H-03 (blue), PYTHIA 8.145 (yellow), SIBYLL 2.1 (green), EPOS 1.99 
(magenta) and DPMJET 3.04 (red). Only the statistical uncertainty of DPMJET 3.04 is 
shown by the error bars as representative of the models. 
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Figure 5: Ratio of the MC spectra divided by the data in each energy bin. The left 
and the right panels show the spectra for the pseudorapidity ranges of r\ > 10.15 and 
8.77 < 77 < 9.46, respectively. The colored plots indicate the results of MC, QGSJET 
h-03 (blue), PYTHIA 8.145 (yellow), SIBYLL 2.1 (green), EPOS 1.99 (magenta) and 
DPMJET 3.04 (red). To describe the size of the errors, the experimental data are also 
shown on the ratio of unity with the statistical uncertainty (error bars) and the total 
uncertainty of data (the black hatches). The luminosity uncertainty was not included. 
The statistical uncertainties of each MC are shown as the error bars of MC. 
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